Post-transient relaxation in graphene after an intense laser pulse 
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Abstract 

High intensity laser pulses were recently shown to induce a population inverted transient state 
in graphene [T. Li et al. Phys. Rev. Lett. 108, 167401 (2012)]. Using a combination of hydro- 
dynamic arguments and a kinetic theory we determine the post-transient state relaxation of hot, 
dense, population inverted electrons towards equilibrium. The cooling rate and charge-imbalance 
relaxation rate are determined from the Boltzmann-equation including electron-phonon scattering. 
We show that the relaxation of the population inversion, driven by inter-band scattering processes, 
is much slower than the relaxation of the electron temperature, which is determined by intra-band 
scattering processes. This insight may be of relevance for the application of graphene as an optical 
gain medium. 
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FIG. 1: Schematic demonstration of the formation of population inverted electronic state right 
after an intense laser pulse, (a) Photoexcited carriers generated by ~10fs pump pulse; (b) The 
leading scattering processes of photo-excited carriers taking place in several femtoseconds: e + e — > 
e + e, h + h — >■ h + h, e + h — >■ e + h (e-electron carriers and h-hole carriers), which quickly establish 
individual thermalization in electron and hole carriers sharing a common electronic temperature T 
due to the electron-hole scattering events; (c) After the internal thermalization, the photoexcited 
carriers form a population inverted hourglass-like electronic state characterized by distinct chemical 
potentials and a common electron temperature T. 

I. INTRODUCTION 

Recently, it was shown that photoinduced femtosecond nonlinear saturation, trans- 
parency and stimulated infrared emission of extremely dense fermions in graphene mono- 
layers emerge pQ. A single laser pulse of 35 fs quasi-instantaneously builds up a broadband, 
inverted Dirac fermion population, where optical gain emerges and manifests itself via a 
negative optical conductivity, see Fig. [I] Increasing the excitation from the linear to the 
highly nonlinear regime, the photoexcited transient state evolves from a hot classical gas to 
a dense quantum fluid. Such high-density population inversion at femtosecond time scales 
has significant implications in advancing graphene based above terahertz speed modulators, 
saturable absorbers, or an ultra-broadband gain medium. These results emerge in a regime 
where the photoexcited carrier density is much larger than the initial background carriers 
and one is no longer in the linear power dependence of transient signals. An important 
open question in this context is the origin of the comparatively stable population inverted 
state despite the rapid thermalization. In Ref. |T] we argued, based on results obtained using 
perturbation theory with respect to the electron-electron Coulomb interaction [2H1], that the 
transient state of dense Dirac fermions is stabilized by the phase-space constraints of the 
relativistic spectrum. The change in the optical conductivity as function of photoexcited 
carriers was then very well described in terms of a nonequilibrium electron distribution 



2 



function. For this distribution function we assumed the quasi-equilibrium form 

characterized by the linear dispersion relation e^x = \vh\k\, where A = + (— ) refers to 
the upper (lower) branch of the graphene spectrum with velocity v. T is the electron 
temperature and n\ refer to the chemical potentials that are allowed to be distinct for the 
upper and lower branch of the spectrum. The population inversion is therefore characterized 
by 5fi = (//+ — //_)/2. Results for the electron temperature and the chemical potentials were 
determined from an analysis of the energy and charge balance of the system [U [S]. 

A natural question to ask is the nature of the post-transient relaxation that gives rise to 
a relaxation of the laser induced population inversion back to equilibrium. In addition to 
the electron-electron Coulomb interaction, the post-transient regime is characterized by the 
coupling of the electronic systems to the lattice, which is expected to lead to a relaxation 
of the electronic energy and population inversion. The investigation of this question is the 
subject of this manuscript. Our theory is a generalization of the approach used in Refs. [Bi- 
ll] to the regime of population inverted initial states. In Refs.jEJ [7] the energy transfer to 
phonons was analyzed as the dominant low-temperature cooling channel of excited electrons 
in graphene without population inversion. Based on the assumption of a rapid thermalization 
of the electron system the change in the electron temperature 

was determined from the electronic heat capacity C and the cooling power Q = ^jp-- 
This approach is justified in a hydrodynamic regime, where the characteristic change of 
the temperature j^y^^ is much smaller than the microscopic relaxation rate r -1 . The 
cooling power Q was then determined from an analysis of the Boltzmann equationj6j Ej. 

Here we generalize this approach and include the corresponding change in the chemical 
potentials p,\ into account, i.e. we analyze the distribution function Eq jl.l| with time de- 
pendent electron temperature and chemical potential: T — > T (t) and fi\ —> fi\(t). This 
enables us to monitor the temporal evolution of the nonequilibrium state that follows the in- 
tense laser pulse and compare the dynamics of the electron heating and population inversion 
8/i (t). In the first part of our theory section we give a summary of the key hydrodynamic 
relations that apply to our system. In a second step we give explicit results for the cooling 
power and imbalance relaxation that are obtained from an analysis of the corresponding 
kinetic equation. 
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II. THEORY 



A. Hydrodynamic considerations 

We analyze the time evolution of the transient state that is characterized by an effective 
electron temperature T (t) and chemical potentials n\ (t) following an intense laser pulse. 
The latter induces electron heating and a population inversion, see Fig. [I] Within a hydro- 
dynamic description, analogous to Refs.[6-8j, we use a quasistatic description. To this end 
we analyze the internal energy £ (T, and the particle numbers N x (T, of the two 

branches of the graphene spectrum: 

dS = C e dT + Y J d£{ Y + ^~ ) d ^ 

X 

dN x = xxdT + K X dn x . (2.1) 

Here we made the assumption that the occupation of the A-th branch of the spectrum 
only depends on its own chemical potential fi x , and not on the chemical potential of the 
other branch. This assumption will be justified later in explicit calculations of the involved 



kinetic processes. In Eq j2.1| we used the heat capacity C, the compressibilities k + and k_ 
of the upper and lower Dirac cone, respectively, as well as the corresponding changes in the 
occupations as function of temperature x+ an d X~ : 



C{T, 

K X (T, /i A ) 



d£ (T, 

dT 
dN x {T,fi x ) 



A 



xM = ?m>. ( , 2) 

The change in energy as function of chemical potential d£ (T, /i_) / dfi X = d£/dfi X \ T , can 
be expressed in terms of xa- The corresponding Maxwell relation is given as dS/dfi x \ T = 
-N x + T X x- 

For quasiparticles with distribution function Eq jl.l| these response functions are given 

as 

j\j r 

C e = 7^ I ( e kA — Hx) fax (1 — fax) , 
T A,A 

K A = ~Y J fa\ (1 - fax) , 

N f 

XX = 7j4 J ( £ kA - fix) fax (1 - fax) , (2.3) 



where = 4 refers to the valley and spin degeneracy of graphene. We use the notation 
(fcyj Ea • • • and L ■ ■ ■ = I (f^p 



Jk a ' ' ' = I 12% Sa ' ' ' an d Jk ' ' ' = / rl^ 2 " ' ' ' ' and se ^ — 1 in the expressions. The 
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time evolution of these quasistatic states is then determined by 

d£ dT ^ dux 

A 

dN x dT dfi X 

-w = Xx m +Kx ^r- (2 - 4) 

Since the evolution should be done under the condition of fixed total charge ^ = ^ A = 
0, we can express the time dependence of the mean chemical potential /Z as || = — — 
Iff' where X = EaXa, k = Ea k a, #X = Ea X Xx, 5k = Y.x Xk x, and /i A = JI + X5/J,. If we 
introduce the charge imbalance A = Ylx AiVx, we finally obtain for the change in energy £ 
and A: 

In order to have explicit expressions for || and ^ we next resort to a kinetic theory. 



B. Analysis of the kinetic equation 

Next we analyze the Boltzmann equation that leads to an energy and imbalance relaxation 
due to the coupling to optical and acoustic phonons. The time evolution of the transient 
state with effective temperature and chemical potentials, caused by the coupling to phonons, 
is then determined by the distribution function /kA (t). The time evolution of the distribution 
function Eq ]l.l| follows from the Boltzmann equation 

<9/kA (t) 

dt 

with collision term: 
2tt 



4a (t) , (2.6) 



'kA — r~ 

n 

2tt 



/ #k k ;a t/kA (1 - /k',A') + (/kA - /k',A') n q>a ] 5 (e k A - W k _ k ',a - ^k',A') 
7k', A' 

/ Skk'" l/k',A' (1 - /kA) - (/kA - /k',A') n q,J ^ (^kA + ^k-k',a - ^k',A') 
Jk'X' 



where g^^j 0, represent the coupling coefficients of the A-band electrons of wavevector k with 
the a- mode phonons to yield a A'-band electrons of k'. Similar to the hydrodynamic analysis 
presented above, we analyze the total energy and the particle numbers of the two branches: 



£ (t) = N d / (e kX - fix (t)) / k A (t) , 
N x (t) = N d jf vx {t). (2.7) 
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From the Boltzmann equation follows 



as 

~dt 
dA 
~dt 



f f 1 

N d / (^kA - A*a) hx -N d / k A- 

7k,A Jk,A 



djix 
dt ' 



N d l A/kA. 

'k,A 



(2- 



The changes in the energy and particle number are now determined by the collision integral 
of electron-phonon scattering. 

We can also make contact between this approach and the hydrodynamic considerations 
presented earlier. The quasi-equilibrium form, Eq jl.l| implies that 



kA 



which yields 



dS 

~di 



dA 

~dt 



dt 



-N d 



kA 



fl 



kAy 



(^kA ~ fJL\) 8T 1 djl X 

' dt T dt 



7^2 

e 



(2.9) 



k.A 



djix 
dt 



kA 



(2.10) 



+N d / /kA (1 - /kA) 



k.A 



(£kA - ^a) dT £ kA - ii\ d/2x 



T' 



dt 



T dt 



N d [ /kA (1 - /kA) A 
Jk,A 



(gkA - A*a) dT 1 djxx 
T 2 dt T dt 



This yields our earlier result Eq 2.4 including the expressions Eq j2.3| for the electronic heat 
capacity, C, the compressibility k x = dN x /dfix and the change of the particle numbers 
with temperature x\ — dNx/dT. This analysis also justifies our earlier assumption that the 
particle number of one branch does not explicitly depend on the other chemical potential. 
Of course, this is a direct consequence of the form Eq ]l.l| of the distribution function. 

We finally obtain a coupled set of equations that determines the time evolution of the 
population inversion 8 /a and electron temperature T: 



C 
f ' 



t\ dT 
K J dt 
8kx\ dT 
k ) dt 



+ [Sx- 



+ K 



8t$x\ d8fi 
88 fi 



K 

8jf 



dt 



N, 



k,A 



£ kA — A 4 A 
T " 

A/kA 



'kA 



(2.11) 



k,A 



In what follows we analyze this set of equations numerically. We present our data as function 
of the typical time scale for r = 47r 2 ft 3 t> 2 / {a\g 1 8ii (0)) with electron phonon coupling con- 
stant g, that characterizes intra-band scattering processes in the regime where the initial pop- 
ulation inversion 8ll (0) is large compared to the phonon frequency. Here, a$ = 1.42 x 10 _10 m 
is the C-C bond length. 
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III. RESULTS 



Our analysis of the time-evolution of the electron temperature T(t) and the population 
inversion S/j, (t) is performed for graphene at the neutrality point. We determine the initial 
electron temperature T (t = 0) and chemical potentials fi\ (t = 0) from the charge and energy 
balance right after the pulse; for details see [H [5]. Our time t = point refers to the 
beginning of the post transient evolution, i.e. about lOOfs after the initial laser pulse that 
caused the population inverted state in the first place. In addition to the value of the Fermi 
velocity, and the equilibrium chemical potentials and temperature, the parameters of the 
theory are the pump-laser frequency hu pump = 1.55 eV and the number of photoexcited 
carriers n ex (t = 0) that is determined from the pump-fluence of the laser. Here we use 
n ex (t = 0) = 4 x 10 13 cm" 2 . We obtain T (0) = 3336.4 K and 5fi (0) = //+ (0) = -//_ (0) = 
0.549761eV. Finally to describe the relaxation due to optical phonons we use hu> ~ 196 
meV as optical phonon energy. For the lattice temperature we assume that it stays constant 
around room temperature Tl = 300K, i.e. we assume that the heat is quickly transferred to 
the bulk of the substrate. This aspect should be more subtle in case of suspended graphene. 
For simplicity we assume that the frequency is momentum and pho non-branch independent. 

The electron-optical phonon coupling constant g 2 = / \ takes a typical value g = 

0.315 eV[10j. Since we present our results as function of t/r, only the overall time scales are 
determined by the value of g. Using 5fi(0) = 0.55 eV we obtain r ~ 10.21ps. The occurrence 
of the various temporal regimes that follow from our theory are not affected by the value of 
9- 

In Fig.[2]|4] we show our results for the changes of the electron temperature T(£), the 
population inversion Sfj,(t), and the density n ex (t) of photoexcited carriers as function of 
time. The density of photoexcited carriers 

n ex = \ (N + - N° + - (JV_ - N )) (3.1) 

results from the calculated chemical potentials. Here N® are the particle densities of the 
two branches of the spectrum before the pulse. 

Three distinct time regimes during the relaxation of a initially hot dense state, with 
/csT(0) > hhj Q > ksTi, are indicated from the numerical result: At short time scale, t/r < 
1/2, the electron temperature drops rapidly, while 5/x rises due to the establishment of 
the sharp Fermi-distribution edge that is associated with cooling. While n ex , decreasing 
slowly and linearly, is not sensitive to the rapid temperature drop and changes similar 
to the population inversion. The fast initial cooling is mainly due to carrier intra-band 
transitions by emitting optical phonons. This regime, dominated by carrier cooling process, 
characterizes the energy relaxation of the system. Next we identify an intermediate regime, 
1/2 < t/r < 10. In this regime we observe a plateau for the evolution of the electron 
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FIG. 2: Time evolution of the electron temperature T(t) for a hot, dense, population-inverted 
electron gas in graphene induced by an intense laser pulse at t < 0. The initial temperature of this 
system, T (0) = 3336.4 K, follows a sharp drop within a short time period t < r/2 mainly due to 
intra-band transitions, then reaches a plateau for a relatively longer time period r/2 < t < lOr 
before further lowering down to the equilibrium temperature. Here r = A^ffiv 2 / (a^g 2 5^{&j) ~ 
10.21 ps. 



temperature, while Sfi (t) and n ex (t) decrease gradually. In this regime, the carrier cooling 
driven by intra-band transitions is less efficient while a relaxation dominated by inter-band 
transitions becomes dominant. This is the regime where the population inversion is being 
destroyed. Particle-hole recombination processes, characterize the population imbalance 
relaxation. Finally, for longer times, t/r > 10 — 20, the relaxation of population inverted 
configurations is essentially completed and inter-band transitions from the upper to the lower 
branch becomes less significant. What is left is a slow cooling by the inefficient intra-band 
transition as ksT(t) <C hw . In this regime, the coupling to acoustic phonons, ignored in 
our treatment, should come into play and eventually become the most dominant relaxation 
process for the terminal relaxation towards equilibrium. 



IV. CONCLUSIONS 



In conclusion, we investigated post-transient state relaxation of hot, dense, population 
inverted electrons in graphene that emerged as the result of an intense laser pulse, as shown in 
Ref . [T| . Using a combination of hydrodynamic arguments and a kinetic theory we determined 
the cooling rate and charge-imbalance relaxation rate. The latter are determined from an 
analysis of the Boltzmann-equation where we included the scattering between electrons and 
optical phonons. We demonstrated that the relaxation of the electron temperature, driven 
by intra-band scattering processes, is much more rapid than the relaxation of the population 
inversion, which is determined by inter-band scattering processes. Thus, the relaxation of 
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FIG. 3: Time evolution of the population inversion 5/j,(t) for a hot, dense, population-inverted 
electron gas in graphene at the charge neutrality point where 5fi(t) = fJ-+(t) = —[i-(t). The initial 
population inversion of this system, dfj, (0) = 0.549761eV, follows a small jump within the short 
time period t < r/2 due to the establishment of the sharp Fermi-distribution edge that is associated 
with cooling, then decreases gradually as particle-hole recombination processes through inter-band 
transitions become dominate for the longer time periodr/2 < t < lOr. Here r ~ 10.21 ps. 
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FIG. 4: Time evolution of the carrier density n ex (t) for a hot, dense, population-inverted electron 
gas in graphene. The photo-excited carrier density takes an initial value at n ex (i = 0) = 4 x 
10 13 cm~ 2 , then drops gradually as a result of the inter-band carrier scattering processes through 
coupling to optical phonons. Here r ~ 10.21 ps. 

the population inversion is significantly slower that the timescales responsible for the energy 
transfers between the hot electron gas and the lattice. This insight may be of relevance for 
the application of graphene as an optical gain medium. 
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